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Abstract. We present spectroscopic analyses of 17 faint 
blue stars in the globular cluster NGC 6752, using optical 
and UV spectrophotometric data and intermediate resolu- 
tion optical spectra. Effective temperatures, surface grav- 
ities, and helium abundances of the stars are determined 
and compared to theoretical predictions. All stars are he- 
lium deficient by factors ranging from 3 to more than 100, 
indicative of gravitational settling of helium. Stars with ef- 
fective temperatures above about 20000 K (sdBs) fit well 
to the evolutionary tracks, whereas the cooler stars show 
lower surface gravities than theoretically expe cted. T his 
agrees with earlier findings by Moehler et al. ( 1995a ) for 
the faint blue stars in M 15. Deriving masses from the at- 
mospheric parameters and the cluster distance leads to a 
mean mass of the sdB stars of 0.50 Mq and a standard 
deviation of about 0.043 dex, which is below the value de- 
rived of the observational errors and therefore consistent 
with a very narrow mass distribution. The mean mass is 
in good agreement with the (0.49 ± 0.02) Mq expected 
if sdB stars are extreme horizontal branch stars with a 
helium core of 0.48 and an extremely thin hydro- 
gen layer. The cooler stars however show a significantly 
lower mean mass of 0.30 Mq with a standard deviation 
of the mean of 0.059 dex in agreement with the findings 
of Moehler et al. ( |1995aD and de Boer et al. ( |1995| ) for 
the BHB stars in M 15 and NGC 6397, respectively. The 
problems presented by the different mass distributions of 
BHB and sdB stars within one cluster are discussed. 
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* Based on observations collected at the European Southern 
Observatory 



1. Introduction 

The horizontal branches (HBs) in colour-magnitude dia- 
grams (CMDs) of galactic globular clusters show a large 
variety in morphology: On the blue side of the RR Lyrae 
gap one may find HBs extending down to very faint mag- 
nitudes at the same blue colour (e.g. M 5, Buonanno et 
al., 1981) while others are short and limited to the region 
0"0 (e.g. NGC 6397, Alcaino & LiUer, 



near My = 0'!'0 (e.g. NGC 6397, Alcaino & LiUer, |198§). 
Also, some of these HBs appear to have gaps (e.g. NGC 
288, Buonanno et al., |1984|; M 15, Buoanno et al. , |1985D 



or, at least, are not populated in a smooth manner (e.g. 
NGC 6752, Buonanno et al. 



1986) 



To distinguish the long and almost vertical extensions 
from the main body of the horizontal branch, the term 
"extended horizontal branch" (EHB) in contrast to "blue 
horizontal branch" (BHB) is used. In the field of the Milky 
Way the EHB stars are identified with the subdwarf B 
(sdB) stars, which have Tes > 20000 K and logg > 5 



1984; Moehler 



(Greenstein & Sargent, 1974; Heber et al. 
et al. [1990| Saffer et al.,^99|). 

Two basic scenarios for the origin of the sdB stars have 
been put forward in the literature: sdBs are assumed to 
have either originated as single stars, which for some rea- 
son have lost more of their hydrogen envelope by the time 
they arrive on the HB than most stars do, or they are 
the produ ct of some kind of binary interaction (see Bai- 
lyn et al. 1992, for a detailed discussion). In the single 
star scenario, the sdB mass is close to 0.5 Mq and de- 
termined by the canonical helium core mass at the core 
helium flash (see Sweigart, 1994). The importance of bi- 



nary evolution is supported by the high binary frequency 
observed amongst the field sdB stars (AUard et al., 1994; 
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Theissen et al., 1994). There are three variants of the bi- 
nary scenario: 

(I) Common envelope evolution can lead to the forma- 
tion of subluminous O and B stars if Roche lobe over- 
flow occurs during the ascend of the first (case B) 
or se c ond (case C) gia nt branch (Iben & Tutukov 



1985| , |l993t Iben |1986| ; Iben & Livio |l993| ) 



WeU 

studied examples are HD 128220 (sdO-|- FIV) and 
AA Dor. For case B evolution the sdB masses are some- 
what smaller than the canonical 0.5 Mq (e.g. AADor, 
M = 0.3 Mq, Kudritzki et al., |1982D , while for case C 
they are slightly larger than the canonical value (e.g. 
HD 128220, M > 0.55 Mq, Howarth & Heber, 



1990) 



(II) Mengel et al. (1976) explored the possibility that sdB 
stars could evolve from binaries in which mass transfer 
is initiated at the tip of the giant branch, i.e. shortly 
before the ignition of the core helium flash, resulting 
in sdB masses very close to the canonical 0.5 Mq. 
sdB stars can also result from the merger of a dou- 
ble helium whi te dwarf system (Iben & Tutukov, 



(III 



1984| ; Iben [1990D with masses ranging from 0.3 Mq to 
0.9 Mq and low mass stars formed preferably. Bailyn 



& Iben ( 198S| ) have explored this possibility for glob- 
ular cluster sdBs and argue that a few tens of sdBs 
can be formed in a typical globular cluster, because 
binary-single star interactions increase the merger rate 
considerably, especially for the central regions of the 
cluster. 

Therefore at least four scenarios for the origin of sdB stars 
are at hand. They can be distinguished observationally 
if the mass distribution of sdB stars can be determined, 
since the four scenarios predict quite different mass distri- 
butions: the single star and Mengel's binary scenario (II) 
predict a sharp mass distribution peaked at the canonical 
mass of 0.5 Mq, while binary scenarios (I) and (III) pre- 
dict broader distributions peaked at masses below 0.5 Mq. 
In order to find out which is the dominant sdB production 
channel in globular clusters, we have started to determine 
spectroscopic masses of blue horizontal br anch s tars in the 
globular clusters M 15 ( Moehl er et al., 1995a paper I), 



NGC 6397 (de Boer et al. Jl995| , paper II), and NGC 6752 



(Heber et al., 1986 and this paper) from atmospheric pa- 



rameters and the known cluster distances. The previous 
analyses of the BHB and EHB stars in NGC 6752, M 15 
and NGC 6397 showed that 

1. Three stars in NGC 6752 are bona fide EHB stars and 
cannot be distinguished from field sdB stars. One very 
blue star is a sdOB star which already evolved from 
the EHB. 

2. The stars below the gap in M 15 are still BHB stars ac- 
cording to their atmospheric parameters Tcff and logg. 

3. The BHB stars in M 15 and NGC 6397 show masses 
that are significantly lower than predicted by canonical 
theory. 



NGC 6752 is up to now the only cluster where the gap seen 
in the CMD really separates sdB stars from BHB stars. 
For a long time it has also been the only cluster where sdBs 
have been securely identified at all. Additional sdB/sdOB 
candidates have been verified recently in M 15 (Durrell 
Harris, |1993| ; Moehler |1995| ) and M 22 (Moehler et al.. 



1995b| ) 



NGC 6752 has a large population of EHB stars 
and is therefore very well suited for our aims. Early 
CMDs (e.g. Cannon, 1981 ) display a pronounced gap 
(between V « 16'°0 and V « 17'°0) separating the EHB 
(V > 17"0) from the BHB (V < 16'^'0). More recent 



CMDs (Buonanno et al. 1986), however, suggested that 



the gap is a region of the HB which is sparsely populated 
rather than devoid of stars. Therefore, we aim to study 
more bona flde EHB stars (below the gap, V < 17™0) to 
determine their masses spectroscopically using the same 
methods as in Heber et al. ( 1986D and paper I. The second 
aim is to clarify the nature of the stars inside the HB gap 
region (i.e. 16™0 < V < 17™0), discovered by Buonanno 
et al. (1986). We also include two stars brighter than the 



HB gap because of their unusually blue colour. 



2. Observations 
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Fig. 1. The colour- magnitude diagram of NGC 6752 as pub- 
lished by Buonanno et al, (1986), Stars analysed in this paper 
are marked by open squares, open stars mark the objects dis- 
cussed by Heber et al, (1986), The star 3-188 of Heber et al, is 
not in the Buonanno data. 
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We selected 17 blue stars from Buonanno et al. (1986) 
for spectroscopic follow-up spectroscopy. These stars are 
marked in Fig. 1 with squares. Open stars in this figure 
mark the stars discussed by Heber et al. ( 19861 ). Nine stars 
lie below the gap in the CMD (V > 17" 0), five stars are 
in the gap region (16™0 < V < 17™0) and three stars are 
above the gap (V < 16™0), two of which have very blue 
colours. All optical spectra were obtained at telescopes 
of the ESO La Silla observatory using focal reducers. For 
calibration purposes we always observed 10 bias frames 
each night and 5-10 flat-fields with a mean exposure level 
of about 10000 ADU each. We obtained low resolution 
spectrophotometric data with large slit widths to analyse 
the flux distribution and medium resolution spectra to 
measure Balmer line profiles and helium line equivalent 
widths. 

The atmospheric extinction was corrected using the 
data of Tiig (1977) and the data for the flux standard 
stars were taken from Hamuy et al. (1992). 



2.1. Observations in 1992 

In 1992 we used EFOSCl at the 3.6m telescope with a see- 
ing between 1.5" and 2". Details of the setup are given in 
Table 1. Since the RCA chip showed some column struc- 
ture we took flat-fields using U and B filters to achieve 
mean count levels between 5 and 10000 counts to allow a 
correction of the column structure. As flux standard star 
we observed Feige 110. During these observations the slit 
angle was kept at its default east-west orientation. 

2.2. Observations in 1993 

In 1993 we observed with EMMI at the NTT, using the 
two-channel mode for the low resolution spectrophotomet- 
ric observations to get a longer wavelength base. Seeing 
values for these observations varied between 0.8"and 1.5". 
Since no EMMI grism allows observations below 3600 A 
(necessary to measure the Balmer jump) we used a grating 
in the blue channel and reduced the dispersion by binning 
along the dispersion axis by a factor of 3. For the inter- 
mediate resolution spectra we used only the blue channel 
and did not bin. the slit was kept at ist default east- 
west orientation {— 0°) except for the spectrophotometric 
observations of B 2395 (75°), B 4009 (75°), and B 4548 
(85°). Positive angles mean anti-clockwise rotation. Since 
we used a narrow slit for the medium resolution obser- 
vations and did not try to correct for atmospheric dis- 
persion there are some slit losses at the blue end of the 
medium resolution spectra. As flux standard stars we used 
LTT 6248, LTT 9491, and Feige 110. 

2.3. Observations in 1995 

In 1995 we used EF0SC2 at the 2.2m MPI/ESO tele- 
scope to re-observe the low-resolution spectrophotometric 
data obtained with EMMI (see Section 3.2). Hence we only 



used the low resolution setup. Here we had seeing values 
between 1" and 1.5". Unfortunately the two nights we 
had were not really photometric. As we had observed the 
standard stars quite frequently we could however check 
the photometric conditions from the response curves. As 
flux standard stars we used EG 274 and Feige 110. Dur- 
ing these observations we did not rotate the slit, but kept 
it at its default east-west orientation. 

2.1 lUE data 

In addition to the observations described above low reso- 
lution lUE SWP spectra are available from the archive for 
the stars B 491, B 916, B 4009, and B 4548. In addition 
to the standard reduction and calibration the correction 
for the new white dwarf scale as described by Bohlin et 
al. ( |1990D and Bohlin ( |1996| ) was performed, we also re- 
calibrated the old lUE data used in Heber et al. ( 1986| ). 
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B 4380 
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X 








-30 


B 4548* 


-307 


+444 
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Table 2. List of observed objects. The star numbers refer to 
Buonanno et al. (1986). Comparing the relative coordinates 
given there to other photometry of NGC 6752, we found that 
X_B = -Y and vice versa. For stars marked with * lUE SWP 
are available. 1 (m) stands for low (medium) resolution spec- 
tra. X marks spectra observed with the respective setup, g) 
marks low-resolution spectra obtained under non-photometric 
conditions. Vhd. is the heliocentric velocity. 



3. Data Reduction 

We describe the reduction of our data sets in some detail 
here since we encountered several problems, esp. with the 
spectrophotometric data, during the reduction and cali- 
bration of our observations. Since these problems on one 
hand affect the quality of our data (and in consequence 
the reliability of our analysis) and were on the other hand 
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Instrument / 
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pixel no. of 
size pixels 
[/im] 


gain read- out 
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fp- / ADTTl [p~1 


EFOSCl 
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15 656 X 1024 


2.1 33 
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1993/07/23-26 NTT 
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["] 




[A/mm] 


[A] 
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5.0 


UV300 
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3200 - 6400 


11 
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EFOSCl/med 


0.75 


B150 
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3800 - 5600 


12 


0.3 


EMMI/DIMD B 


5.0 


#4 


72 


3400 - 5300 


7 


0.8 


EMMI/BLMD B 


1.0 


#4 


72 


3200 - 5100 


15 


0.3 


EMMI/DIMD R 


5.0 


# 13 


224 


4000 - 8500 


14 


0.5 


EFOSC2 


5.0 


# 1 


442 


1000 - 9400 (nominal) 
3400 - 9400 (used) 


16 


1.7 



Table 1. Observation and reduction parameters 



produced by instruments much used by observers at ESO 
we think that a detailed report is adequate. 

3.1. EFOSCl and EMMI data 

We always averaged bias and flat-field frames for each 
night separately, since there were small variations in the 
bias frames (up to a few counts/pixel) and also a slight 
variation in the fringe patterns of the flat fields from one 
night to the next (below 5%). The dark currents were de- 
termined from several long dark frames and turned out 
to be negligible in all cases. To correct the electronic off- 
set we scaled the bias frames by the mean overscan value 
of the science frame. To take the spectral signature out 
of the averaged flat field frames we fitted a 5**^ (7'^) or- 
der polynomial to the blue (red) channel EMMI data. The 
EFOSCl flat fields could not be fitted well enough by poly- 
nomials. We therefore averaged the mean flat-field of each 
night along the spatial axis, smoothed it heavily to erase 
all small scale structure and divided the original mean 
flat by the smoothed average. For the EFOSCl data we 
also had to correct column structure. This correction was 
performed as described in paper I. 

For the wavelength calibration we fitted 3'''^ order 
polynomials to the dispersion relations. We rebinned the 
frames two-dimensionally to constant wavelength steps. 
Before the sky fit the frames were smoothed along the 
spatial axis to erase cosmics in the background. To deter- 
mine the sky background we had to find regions without 
any stellar spectra, which were sometimes not close to the 
place of the object's spectrum. Nevertheless the flat field 
correction and wavelength calibration turned out to be 
good enough that a linear (EFOSCl) or constant (EMMI) 
fit to the spatial distribution of the sky light allowed to 
subtract the sky background at the object's position with 



sufficient accuracy. This means in our case that we do not 
see any absorption lines caused by the predominantly red 
stars of the clusters or by moon light and that the night 
sky emission lines in the red part could be corrected down 
to a level of a few percent. The fitted sky background was 
then subtracted from the unsmoothed frame and the spec- 
tra were extracted using Home's algorithm (Horne, 1986 ) 
as implemented in MIDAS. 

Finally the spectra were corrected for atmospheric ex- 
tinction and response curves were derived from the spectra 
of the flux standard stars. Normally the response curves 
were fitted by a 3'''^ order spline, except the red EMMI 
data, for which we used a 6'^ order polynomial. We took 
special care to correctly fit the response curves for the blue 
low resolution data in the region of the Balmer jump, since 
we will use this feature for the determination of the effec- 
tive temperatures of the stars. 

Comparing the EFOSCl and EMMI low resolution 
data for the stars observed during both runs showed that 
the EFOSCl data have significantly less UV flux blueward 
of the Balmer jump than the EMMI data. Since both sets 
of spectra could be fitted with Kurucz (1992) models (but 
yielded temperatures different by up to 10000 K) we could 
not decide which data set was the correct one and whether 
one or even both data sets were hampered by instrumental 
effects that were not removed by the described reduction. 
We discussed these problems extensively with people at 
ESO and finally decided to repeat the spectrophotomct- 
ric observations with a third instrument (which decision 
resulted in the 1995 observations with EF0SC2). 

The medium resolution spectra observed in 1992 and 
1993 are plotted in Fig. 2 and 3. Besides the Balmer lines, 
He I lines can easily be identified in the spectra of several of 
the programme stars. In addition Hen, 4686 A, is visible 
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Fig. 2. Normalized medium resolution EFOSCl spectra. The part shortward of 3900 A was normalized by taking the highest 
flux point as continuum value. The spectra are sorted along decreasing Balmer line depths. Small tickmarks mark steps of 20%. 



in B 852, which therefore has to be classified as spectral 
typ e sdO B according to the scheme of Baschek & Nor- 
ris ( 1975 ). It is also worthwhile to note that weak Mgii, 
4481 A, absorption is detected in B 1509 and B 3975 

We also used the medium resolution data to derive 
radial velocities, which are listed in Table 2 (corrected to 
heliocentric system). The error of the velocities is about 
30 km/sec. Within error limits the radial velocities agree 
with the cluster velocity (—32.1 km/s, Pryor & Meylan, 



19931) 



3.2. EF0SC2 data 

For the EF0SC2 data we averaged the bias and flat field 
frames over both nights of the run, as they showed no 
deviations above the 1% level. The bias correction was 
again performed by scaling the bias with the mean over- 
scan of the frame and we did not correct for dark current. 
The normalization of the flat-field was somewhat difficult: 
We fitted a 7'^ order polynomial to the flux distribution 
of the averaged flat fields, which did not fit well enough 



the UV/blue part. We therefore averaged and smoothed 
the flat as described for the EFOSCl data and merged 
the UV/blue part of the smoothed averaged flat and the 
red part of the fit at a position where the two overlapped 
and had a similar slope. This merged frame was smoothed 
again before correcting the flat. 

For the wavelength calibration we again used a 3'^^ 
order polynomial to fit the dispersion relation. However, 
while the wavelength calibration frames showed maximum 
deviations of 4 A when rebinned using this dispersion rela- 
tion, the science frames exhibited systematic shifts to the 
red in the order of 30 - 40 A. The only difference between 
the two types of frames (beside the illumination) was the 
slit width, which was 1" for the calibration frames and 5" 
for the science frames. A possible explanation of this effect 
could lie with the fact that the polynomial fits normally 
used for wavelength calibrations are only an approxima- 
tion to the true shape of the dispersion relation, which 
is described correctly by trigonometric functions (Strocke 
1967 ; Bahner 1967 ). Tests performed by Rosa & Hopp 
(1995, priv. comm.) showed, however, that fits using such 




Fig. 3. Normalized medium resolution EMMI spectra. The part shortward of 3900 A was normalized by taking the highest flux 
point as continuum value. The spectra are sorted along decreasing Balmer line depths. Small tickmarks mark steps of 20% . 
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trigonometric functions are extremely unstable unless a 
very good guess of the true relation is already available. 
We therefore decided to correct the offsets by applying a 
mean shift of 36 A to the blue. 

Sky subtraction, correction of atmospheric extinction, 
and flux calibration were performed in the same way as 
for the EMMI and EFOSCl data, using a constant for the 
spatial profile of the sky background and a spline fit for 
the response curve. As mentioned above we could judge 
the photometric quality of the observations from the re- 
sponse curves. We used only those curves that yielded the 
maximum response since those should result from observa- 
tions taken without any obscurations. Several object spec- 
tra calibrated with these response curves showed an offset 
from the B and V magnitudes taken from Buonnano et al. 
( |l986|) , being about 0.2 mag too faint. Although there was 
hardly any wavelength dependency of this offset visible we 
could not be sure that the absorption was truly grey. We 
therefore decided not to use these data for determinations 
of the effective temperatures because we could not rely on 
the continuum slope. They are marked in Table 2 by (g). 



3.3. Comparison of the optical spectrophotometric data 

Comparing the data for the three runs shows that the 
EMMI and the EF0SC2 data agree rather well when a cor- 
rection for the much lower resolution of the EF0SC2 data 
is applied (necessary to allow comparison of the Balmer 
jump region). This leads to the assumption that there 
is some problem with the UV response of the EFOSCl 
that is not corrected by flux calibration. An additional 
reason could lie with the extraordinarily high extinction 
observed at La Silla in the years 1991 and 1992 (Burki 
1995| ). We therefore ignore all low resolution data 



et al. 

obtained with EFOSCl for the further analysis. As a con- 
sequence, 5 (B 210, B 852, B 1509, B 1628, B 4380) out of 
17 programme stars lack reliable spectrophotometric data. 
For an additional star (B 3915) we do not have any low 
resolution data at all. 



4. Atmospheric Parameters 

Effective temperatures, gravities, and helium abundances 
are derived from line blanketed LTE model atmospheres. 
We used the Kurucz ( 1992 ) grid and models calculated 
with an updated version of the code of Heber (1983) 
using ATLAS6 opacity distribution functions (ODF), 
log(Z/ Zq) — —2.0 and helium abundance of 0.1 solar, ap- 
propriate for blue horizontal branch stars . 

Kurucz (|1992|) ATLAS9 models (logg =5.0, solar He 
abundance, log(Z/ Zq) = — 1.5 closest to the cluster's 
metallicity, [Fe/H] = —1.54, Djorgovski, 1993) were used 
to fit the energy distributions (Sect. 4.1) because of their 
superior spectral resolution in the region of the Balmcr 
jump. To correct for interstellar reddening we applied the 



reddening law of Savage & Mathis (1979) and used an 
Eb-v of 0.05 mag. 

Heber's models were used to analyse the hydrogen 
and helium line spectrum (Sect. 4.2 and 4.3) in order to 
account for the appropriate gravities and helium abun- 
dances. As these models were computed for [M/H] = —2.0 
(instead of —1.5) and sub-solar helium abundance we fit- 
ted the Kurucz hydrogen hue profiles (Teff = 22000 K, 
logg = 5.0) by Heber models of log(Z/ZQ) = -2.0. The 
use of different generations of ODFs (ATLAS9 versus AT- 
LAS6) and different helium abundances turns out to be 
neghgibly small, both for Teff and logg. Hence, the fol- 
lowing procedure is applied: The Balmer lines Hp, H-^, 
and Hs are fitted using Heber's models. At the final pa- 
rameters (Teff, logg) curves of growths for the Hei hues 
4026, 4121, 4388, 4471 and 4921 A are constructed to de- 
termine the helium abundances from equivalent widths of 
these lines. 

In order to allow a consistent treatment of the results 
of Heber et al. (1986), which were obtained with solar 
metallicity models, we checked the effect of the metallic- 
ity on Teff and logg by fitting solar metallicity models 
with the [M/H] = —2.0 models we use now. We found out 
that the solar metallicity models yield effective tempera- 
tures about 1000 K lower than the metal poor models, but 
essentially the same logg values. We therefore increased 
the Teff values of Heber et al. (1986) by an offset of 1000 K 
to account for the lower cluster metallicity and kept the 
logg values. 

In the case of B 852 LTE models were found to be in- 



appropriate due to its high Teff (see Napiwotzki, |1996|). 
We used his NLTE models which were calculated with 
the code of Werner & Dreizler ( 1996[ ) and do not include 
metal line blanketing. Therefore the Toff scale slightly dif- 
fers from the metal line blanketed scale of Kurucz and 
Heber models. No attempt was made for adjustment of 
Teff scales. 

4.1- Toff derived from energy distributions 

4.1.1. lUE data and the results of Cacciari et al. (1995) 



For the four stars which have lUE spectra we tried to 
simultaneously fit the continuum data (lUE spectra, BV 
photometry, continuum of the low resolution optical spec- 
tra) and the Balmer jump (cf. Fig 4). In two cases the 
Balmer jump suggested temperatures different from those 
derived from the continuum data (B 491 and B 916, cf. 
Table 3). 

The same data have already been analysed by Cacciari 
et al. (1995) by calibrating UV- (1300A -I8OOA) and UV- 
visual (I8OOA - V) colors in terms of Toff and logg using 
Kurucz model atmospheres. In all cases our effective tem- 
peratures are lower than the ones given by Cacciari et al. 
(1995). This can be traced back to the use of different lUE 
flux calibrations. We used the most recent calibration of 
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2000 4000 6000 




Fig. 4. Fits to low resolution data of B 491, B 916, and B 2395. The models are ATLAS9 (Kurucz, 1992) models of metallicity 
— 1.5, all for logg = 5.0. The asterisks mark the B and V fluxes derived from Buonanno et al. (1986) using the conversion 
factors given by Heber et al. (1984). The fluxes as well as the observed spectra are dereddened using the reddening law of Savage 
& Mathis (1979). 



Bohlin (1996) which gives a flux at I8OOA lower by 10% 
than one of the Cacciari et al. (1995) and, hence, chang- 
ing the (18- V) colour by 10% , which translates in a Tcfi 
decrease exactly as found when comparing our results to 
those of Cacciari et al. 

In view of this flux calibration problem, we also went 
back to the lUE data for the stars analysed by Heber 
et al. (1986) and fltted the recalibrated lUE spectra and 
the optical photometric data as described above. The new 
temperatures were always within 1000 K of the old values 
and no systematic trend showed up. We therefore decided 
to keep the old values and only adjust them for metallicity 
effects (see above). 

4.1.2. Optical data only 

The simultaneous fitting of the spectrophotometric data 
(Balmer jump and continuum) and the B and V mag- 
nitudes worked well (cf. Fig. 4) in all cases except four: 
B 1628, B 2395, B 3655, and B 2126. B 1628 shows a 
strong red excess when compared to the B and V mag- 



nitudes and to the model fitting of the Balmer jump; the 
observed spectrum is also significantly brighter than the 
B and V magnitudes. The stars B 2395 (cf. Fig. 4) and 
B 3655 both show a moderate red excess when compared 
to the model fitting the Balmer jump, which also repro- 
duces the B and V magnitudes. In the case of B 2126 
the optical spectrum agrees well with BV photometry, but 
both show a small red excess when compared to the model 
fitting the Balmer jump. 

Cautioned by these findings we searched all stars for 
red neighbours. For this purpose we used the photometry 
of NGC 6752 by Buonanno et al. ( 1986] ), and extracted for 
each of our targets all stars within a radius of 20" (down 
to the hmiting magnitude of about B w 18™5). We then 
estimated the straylight provided by these stars assum- 
ing that the seeing at larger distances is best described 
by a Lorcntz profile and scaling the intensities with the 
V fluxes. As seeing value we used 1.5"as "intrinsic see- 
ing" for all observations (which overestimates the seeing 
near the meridian) and took into account the elongation 
caused by atmospheric dispersion. Under these assump- 



S. Moehler ' et al.: Hot HB stars in globular clusters - physical parameters and consequences for theory 



9 



tions we get stray light levels of more than 3% in V for the 
following objects: B 617, B 1288, B 1628, B 2162, B 2395, 
and B 3655. As expected we found the largest values for 
B 1628, B 2395, and B 3655, which were observed at rel- 
atively large zenith distances and have very bright neigh- 
bours. For those objects where the neighbour did not lie 
in the slit the observed stray light levels are lower than 
the calculated ones by a factor of 1.5 to 2. The calculated 
excesses for B 1288 and B 617 are about 7% and cannot be 
measured in our spectrophotometric data. Hence we can 
explain the above mentioned red excesses as being caused 
by stray light from nearby stars (except B 2126). At the 
same time we can derive an upper limit of 7% for the 
accuracy of our spectrophotometric calibration (at V) . 

Since the Balmer lines (H^, H^, H^) are diagnostic 
tool we also checked the stray light level in the medium 
resolved spectra. We used the same assumptions as above, 
but scaled the intensities with the B magnitudes. This re- 
sults in straylight levels of at most 5% for B 617, B 1288, 
B 1509, B 1628, B 2162, B 2395, and B 3655. In order to 
estimate the effect of any additional continuum on the 
physical parameters we subtracted from the spectrum of 
B 2395 a constant offset of 5%, renormalized the spec- 
trum and derived ToS and logg again from only fitting 
the Balmer line profiles. It turned out that the best fit 
yielded an effective temperature of 1000 K less than be- 
fore and the surface gravity remained the same. 

4-2. Balmer line profile fits 

To check the validity of the effective temperatures ob- 
tained from the low resolution data we also determined 
Tcff and log g simultaneously by fitting the line profiles of 
H/3, H^, and H5. 

In most cases where low and intermediate resolution 
spectra are available the effective temperatures derived 
from both agree within the error limits (cf. Fig. 5 and 
Table 3). In one case (B 617) however, we find an effec- 
tive temperature from the low resolution data of about 
33000 K, while the Hne profiles yield 25000 K. Since this 
star is in colour and brightness comparable to B 1288 and 
B 2162, we assume that the higher temperature is the valid 
one. The systematic offset between Teff derived from low 
resolution data and from the Balmer line profiles reported 
in paper I could not be found with these data. 

In order to get an estimate for the internal errors of 
the physical parameters derived from fitting the Balmer 
lines we used the independent line fitting code of R. Saffer 
(see Saffer et al., 1994) to verify our results. It turned out 
that we got the same values for spectra with good S/N. 
For the spectra of 1992, which are considerably more noisy 
than the ones from 1993, we get systematically lower tem- 
peratures with Saffer's routine than from our line profile 
fits (using the same set of models). Since we cannot rec- 
oncile the spectrophotometric data with the temperatures 
found by Saffer's routine we decided to keep our values. 
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Fig. 5. Fits to medium resolution data of B 491, B 916, and 
B 2395. The models are those of Heber (1983) for those pa- 
rameters which best fitted the Balmer lines (Toff is already 
transformed to the Kurucz scale): 
B 491: Toff = 28000 K, logg = 5.2 
B 916: Toff = 30000 K, logg = 5.6 
B 2395: Teff = 22000 K, logg = 5.0 



The mean internal error in Teff provided by Saffer's fits 
are 1370 K (1992) and 460 K (1993); the mean internal 
error in logg is 0.17 dex (1992) and 0.06 dex (1993). 



4.3. The case of B 852 

B 852 is the only star in our sample which shows Hen, 
4686 A. Comparing its strength to that of Hei, 4471 A 
suggests a temperature exceeding 35000 K. At such high 
temperature deviations from LTE become important even 
for the relatively large gravities of hot subdwarf stars as 
demonstrated by Napiwotzki (1996). Therefore we used 
his NLTE models to derive Tcff, logg, and helium abun- 
dance from H^, H.^, and H^, He ft, 4686 A and Hci, 4471 A, 
the latter line ratio being a very sensitive Teff indica- 
tor. The best fit is shown in Fig. 6. The depth of the 
observed Balmer line profile cores cannot be reproduced 
by the NLTE models, which we attribute to the lack of 
metal line blanketing. The helium lines are not affected 
by this effect and are reproduced very well. We therefore 
use the best fit values of Toff — 40000 K, logg = 5.2 and 
log(nHc/nH) = —2.0 as atmospheric parameters for B 852. 
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Star 


V 


B-V 




TefF , opt . 


Teff ,lines 


Teff fin 


logg 


M 


Mv 




[mag] 


[mag] 


[K] 


[K] 


[K] 


[K] 




[Mq] 


[mag] 


Blue HB stars 


B 577^ 


14.85 


—0.04 








13700 


3.9 


0.43 


1.61 


B 1509 


15.52 


-0.06 






17000 


17000 


4.1 


0.26 


2.28 


B 2454^ 


14.27 


-0.06 








10700 


3.5 


0.45 


1.03 


B 4104 


15.20 


+0.00 








17000 


4.0 


0.27 


1.96 


o 4 / iy 


it). DO 


— U.UO 








ioUUU 


A A 

4.U 


A OA 

U.2U 


z.4i 








Stars in the 


gap region 










B 1628 


16.30 


-0.17 






21000 


21000 


4.7 


0.35 


3.06 


B 2395 


16.73 


-0.24 




21000 


22000 


21500 


5.0 


0.46 


3.49 


B 3655 


16.40 


-0.22 




24000 


22000 


23000 


5.1 


0.70 


3.16 


B 3975 


16.27 


-0.22 




21000 


21000 


21000 


4.8 


0.46 


3.03 


B 4548 


16.60 


-0.14 


22000 


23000 


23000 


22500 


5.2 


0.75 


3.36 


EHB stars 


B 210 


17.04 


-0.31 






27000 


27000 


5.6 


0.94 


3.80 


B 331^ 


17.12 


-0.22 








26000 


5.6 


0.93 


3.88 


B 491 


17.45 


-0.31 


28000 


30000 


28000 


28500 


5.3 


0.29 


4.21 


B 617^ 


17.84 


-0.33 




33000 


26000 


33000 


6.1 


0.98 


4.60 


B 763^ 


17.13 


-0.24 








27000 


5.5 


0.68 


3.89 


B 916 


17.61 


-0.26 


27000 


30000 


30000 


28500 


5.4 


0.32 


4.37 


B 1288 


17.90 


-0.31 




29000 


27000 


28000 


5.5 


0.31 


4.66 


B 2126 


17.28 


-0.20 




31000 


28000 


29500 


5.7 


0.79 


4.04 


B 2162 


17.88 


-0.27 




35000 


32000 


33500 


5.9 


0.57 


4.64 


B 3915 


17.16 


-0.23 






31000 


31000 


5.5 


0.47 


3.92 


B 4009 


17.44 


-0.30 


33000 


29000 


31000 


31500 


5.7 


0.61 


4.20 


3-118^ 


17.76 


-0.26 








24500 


5.2 


0.23 


4.52 


Post-EHB stars 


B 852 


15.91 


-0.28 








39000 


5.2 


0.53 


2.67 


B 1754^ 


15.99 


-0.24 








40000 


5.0 


0.29 


2.75 


B 4380 


15.93 


-0.14 






32000 


32000 


5.3 


0.93 


2.69 



Table 3. Physical parameters of the observed stars. The star numbers and V, B-V data are taken from Buonanno et al. 
(1986), except for 3-118, whose data were obtained from Caloi et al. (1986). 

^: Data from Heber et al. (1986), adjusted to the Kurucz temperature scale. Note that the mass value for B 1754 given by 
Cacciari et al. (1995) results from a incorrect bolometric correction (Cacciari, priv. comm.) 

^: Since B 617 lies close to B 1288 and B 2162 in the CMD we believe that the higher temperature is the correct one (see text). 
It was not taken into account for the mass distribution. 



4-. 4- Final parameters 

As effective temperature we finally used the mean value of 
all available determinations, assessing the lUE tempera- 
tures double weight because of the superior sensitivity to 
Toff of UV data. Besides the internal errors quoted above 
we also have to account for systematic errors caused by 
the imperfections of models and data. These are estimated 
from the differences between various Tcff determinations 
for the same object to be about 7% in Tes (cf. Table 3). 

For the given final temperature we take that value of 
log g as surface gravity that yields the smallest errors for 
the Balmer line fits at this temperature. An error in Tcff of 
7% translates into an error in log g, which we estimated for 
the three groups of stars mentioned in Table 3: 0.1 dex for 
EHB stars (<Teff> = 15000 K); 0.15 dex for stars in the 
gap region (<Toff> = 22000 K) and 0.2 dex for EHB stars 
(< Teff> — 29000 K). Added to these errors in logg are the 
internal errors given by Saffer's fitting routine (see Sect. 



4.2). For the data given by Heber et al. ( 1986 ) we used 
the mean of the errors for the 1992 and 1993 data. This 
results in the following errors in logg: BHB stars (1992: 
0.196; 1993: 0.115); stars in the gap region (1992: 0.225; 
1993: 0.160); EHB stars (1992: 0.261; 1993: 0.208). Ab- 
solute visual magnitudes were derived from the apparent 
magnitude and the distance modulus ((m— M)v — 13™24, 
Djorgovski 1993). Note that th reddening free distance 
modulus of 13™12 of Djorgovski agrees within the error 
limits with the new distance modulus (m— M)v,o = 13™05 
derived by Rcnzini et al. (1996) by fitting the white dwarf 
sequence. The results are listed in Table 3. 

4-. 5. Helium abundances 

From the medium resolution spectra we could also de- 
termine the helium abundances of the stars (in case he- 
lium lines are seen) or at least upper limits. The mea- 
sured equivalent widths of the Hei lines are given in Ta- 
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edge of the HB in this part. They are probably in a Post- 
EHB stage of evolution (cf. Fig. 7), as already argued by 
Heber et al. ( |1986D in the case of B 1754, and will evolve 
directly towards the white dwarf cooling sequence. The 
distribution of stars in the (Tes, My) diagram (Fig. 8) is 
consistent with that in the {TcS, logg)-plane (Fig. 7) indi- 
cating that the stars above 20000 K (except the three men- 
tioned above) are bona fide EHB stars. The twelve stars 
below the gap (V > 17™0) have a mean absolute mag- 
nitude of My = 4™2 ± 0"^3 in excellent agreement with 
the mean absolute magnitu de of s outh galactic pole sdBs, 

19861 ). The five stars in the 



M 



4-2 ± 0-7 (Heber, 



gap region (16 "0 < V < 17 - O) have Teff > 20000 K, logg 
near 5 and subsolar helium abundances. Therefore they 
show the same characteristics as the field sdB stars and 
have to be classified as sdBs too. However, they are some- 
what more luminous (My = 3™2 ± 0™2) than the stars 
below the gap in NGC 6752 and the field sdBs. 

The fact that the stars between 10000 K and 20000 K 
lie systematically above the HB region in the (Tea, logg) 
diagram has already been noted in papers I and II and 



also by Crocker et al. (1988) 



Fig. 6. Fit to medium resolution data of the sdOB star B 852. 
Overplotted is a Non-LTE model computed by R. Napiwotzki 
(priv. comm.) for a metal-free stellar atmosphere with a 0.1 
solar helium abundance. 



ble 4. Only one star (B 852) shows also Hen 4686 A. Using 
the physical parameters listed in Table 3 we then derive 
helium abundances. The spectra observed in 1992 were 
rather noisy and therefore equivalent widths of less than 
0.35 A cannot be measured. The data taken in 1993 are of 
much better quality and allow equivalent widths as low as 
0.15 A to be measured. We take these values (0.35 A resp. 
0.15 A) as upper limits if we cannot detect a helium line 
in a spectrum and derive an upper limit for the helium 
abundance from the absence of Hei, 4471 A, the helium 
line predicted to be strongest in the observed wavelength 
range. This upper limit is determined for each object in- 
dividually. All stars are helium deficient with respect to 
the sun by factors ranging from 3 to more than 100 and 
no trends with Tcff or logg are apparent, si milarly to the 



TAHB/ 



field sdBs (as demonstrated by Schulz et al. , |199l| ) 



4-6. Evolutionary status 

The results are compared in Fig. 7 {TcS, logg) and 8 (Toff, 
My) to the predictions of stellar evolution theory (Dorman 
et al. 1993). 

The bulk of the stars hotter than 20000 K lies within 
the theoretically predicted region (within the error bars), 
except three (B 852, B 1754, and B 4380). The latter three 
stars have about the same magnitude (V — 15™91-15™99) 
and lie above the gap in the CMD, albeit far to the blue 



Heber et al. (IG 
this paper 



— 5 s 



ZAIIB 



4.6 4.4 4.2 

log T,„ 



Fig. 7. The resulting physical parameters listed in Table 3 
compared to theoretical predictions. The solid lines are the 
Zero Age HB and the Terminal Age (core helium exhaus- 
tion) HB for [Fe/H] = -1.48 of Dorman et al. (1993). The 
long dashed line marks the gap region seen in the CMD by 
Buonanno et al. (1986). The short dashed line gives the posi- 
tion of the helium main sequence (Paczynski 1971). The long 
dashed-short dahed lines give post- EHB evolutionary tracks by 
Dorman et al. (1993), labeled with the total mass of the EHB 
star. 



5. Masses 

Knowing Tcff, logg, and the distances of the stars we can 
derive the masses as described in paper I: 
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Name 


year 


4UzD A 
rXi 


A 1 on A 
4iz0 A 

rXi 


A'\ A A A 

4144 A 

rXi 
lAJ 


/I ooo A 
4388 A 
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44 r 1 A 

rXi 
[AJ 


A no 1 A 
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rXi 
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log UHe 


nHe,0/nHe 


B 210 


1992 


0.75 








1.06 




—2.0 


10 


B 491 


1992 














< —3.2 


>160 


B 617 


1992 














< —2.4 


25 


B 852^ 


1993 










0.26 




—2.0 


10 


B 916 


1993 


0.90 


0.07 


0.30 




1.40 


0.40 


— 1.8 


4 


B 1288 


1992 










1.57 


1.00 


— 1.5 


3 


B 1509 


1992 


0.34: 
















2 


1993 


0.46 




0.12* 




0.45 


0.18* 


—2.2 


15 


B 1628 


1993 


0.28 


0.14 


0.26 


0.22* 


0.38 


0.31 


—2.4 


25 


B 2126 


1992 














< —3.0 


>100 


B 2162 


1993 


0.49 








0.80 


0.45 


— 1.8 


6 


B 2395 


1993 


0.79 






0.73 


0.90 


0.35 


— 1.9 


7 


B 3655 


1992 


0.80 








0.74 


0.95 


— 1.8 
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B 3915 


1993 
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>100 
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n 1 
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10 


B 4009 


1993 














< -3.1 


> 120 


B 4380 


1992 










0.38: 




-2.3 


20 


B 4548 


1993 


0.87 


0.36: 


0.97: 




0.95 




-1.9 


8 



Table 4. The equivalent widths of the He i lines as measured from the medium resolution spectra and the corresponding helium 

abundances (in particle numbers). 

^ B 852 shows also Hen at 4686 A (Wa = 0.8 A) 

^ B 1509 (1993) shows also Mgii 4481 AA (Wa = 0.23 A) 

^ B 3975 shows also Mgii 4481 A (Wa = 0.26 A) and Hei 4169 A(Wa = 0.13: A) 
: equivalent width has a large error due to a rather noisy spectrum 
* equivalent width is at the level of the measuring error 
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Fig. 8. The absolute V magnitudes and effective temperatures 
listed in Table 3 compared to theoretical tracks by Dorman et 
al. (1993, details see Fig. 7). The long dashed line marks the 
gap region seen in the CMD by Buonanno et al. (1986). 
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Fig. 9. The resulting logarithmic masses listed in Table 3 vs. 
log Toff. The objects that were not used to derive the mean sdB 
mass (B 617, B 852, B 1754, B 4380) are marked by circles. 
Also plotted is the Zero Age HB for [Fe/H] = -1.48 of Dorman 
et al. (1993) The long dashed line marks the gap region seen 
in the CMD by Buonanno et al. (1986). 



log (M/M©) = const. + log g + ^A-{(ra-M)y-V^Vth) 
{Yth denotates the theor etical brightness at the stellar sur- 
face as given by Kurucz 1992| .) 

The results arc listed in Table 3 and plotted in Fig. 9. 
In addition to the error in log g errors in the absolute mag- 



also enter the final error in log M. We assume an error 
in the absolute brightness of 0'"l for stars with V < 17™0 
and 0™13 for fainter stars. The error in the theoretical V 
brightness is dominated by the errors in Toff and amounts 



nitude and the theoretical brightness at the stellar surface to 0™12 (BHB stars); 0™13 (stars in the gap region); and 
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0^14 (EHB stars). Alltogether this leads to the follow- 
ing errors in log M: BHB stars (1992: 0.210; 1993: 0.137); 
stars in the gap region (1992: 0.238; 1993: 0.178); EHB 
stars (1992: 0.275; 1993: 0.225). To determine the m ean 
mass of the sdBs (this paper and Heber et al., p_986 ) we 
took all stars with T^s > 20000 K, excluding B 617 (due 
to the strange offset between temperatures from contin- 
uum and from line profiles) and the post-EHB stars B 852, 
B 1754, B 4380, resulting in a total of 16 stars. We then 
calculated the weighted mean of the logarithmic masses 
for these stars (the weights being derived from the in- 
verse errors). The mean logarithmic mass for these stars 
then is -0.300 (= 0.500 M©). The logarithmic standard 
deviation is 0.043 dex and the expected mean logarith- 
mic error as derived from the observational errors is 0.054 
dex. The mean mass therefore agrees extremely well with 
the value of 0.488 Mq predicted by canonical HB the- 
ory for these stars (Dorman et al., 1993) and the stan- 
dard deviation is less than expected. If we omit the stars 
observed in 1992 (due to their higher errors) we get a 
mean logarithmic mass for the remaining 12 stars of -0.327 
(= 0.47 M0), which is somewhat lower than the result 
above but still in very good agreement with theoretical 
predictions. There is no significant difference between the 
mean mass for the eleven stars below the gap (<log M> 
= -0.306; <M> = 0.494 Mq) and the five stars inside 
the gap region (<log M> = -0.289; <M> = 0.514 Mq); 
again pointing towards their nature being identical. The 
five stars above the gap (V < 16™0, Toff < 20000 K) have 
a mean logarithmic mass of —0.518 {— 0.303 Mq) with 
a logarithmic standard deviation of 0.059 (compared to 
an expected error of 0.073). Canonical theory would pre- 
dict for these stars a mean logarithmic mass of —0.256 
(= 0.554 Mq) with a scatter of 0.018 dex. 

6. Conclusions 



tude than the average field sdBs. Their mean mass is 
also consistent with the canonical mass. Their position 
in the (Tea, logg)-plane agrees less well with HB evo- 
lutionary calculations of Dorman et al. (1993) than the 
fainter EHB stars do. 

(v) The stars brighter than the HB gap have significantly 
lower gravities than predicted by theory. Their masses 
are also significantly lower than predicted by evolu- 
tionary calculations. 

(vi) The standard deviations of the mean logarithmic 
masses for all groups of stars are lower than expected 
from the observational errors. 



(i) 
(ii) 



Combining our results with those of Heber et al. (1986) 
we have carried out spectroscopic analyses of 25 blue stars 
in NGC 6752. We studied stars both brighter and fainter 
then the gap region (16™0 < V < 17™0) and also five stars 
within this gap region. The results can be summarized as 
follows: 

(i) All stars are helium deficient and no trends of the he- 
lium abundance with atmospheric parameters become 
apparent. 

(ii) All stars below the gap (except three) lie on or close 
to the Horizontal Branch. The latter three stars have 
already evolved off the EHB towards the white dwarf 
cooling sequence. 

(iii) The stars fainter than the gap (V > 17'"0) are very 
similar to the field population of sdB stars. Their mean 
mass is consistent with the canonical value of half a 
solar mass. 

(iv) The stars inside the gap region (16™0 < V < 17™0) are 
sdB stars too, although more luminous by one magni- 



(iii; 



From these results we conclude: 

The helium deficiency found in all stars is indicative 
for diffusion, i.e. gravitational settling of helium as is 
the case for the field sdB and BHB stars. 
The increased sample size allowed a more meaningfuU 
determination of the mean mass of the sdB stars than 
pevious work did and it turned out to be in good agree- 
ment with the value predicted by classical EHB theory. 
In addition, the standard deviation is smaller than the 
mean error derived from the observational errors. We 
take this as a verification that sdB stars indeed show 
a very narrow mass distribution, which is in full agree- 
ment with their status as extreme HB stars with a 
very thin hydrogen layer, produced by the same pro- 
cesses that produce BHB stars. This finding, however, 
does not explain the gap that separates BHB from sdB 
stars in NGC 6752 (and also in M 15, Durrell & Harris 
1993 ). The mean mass and the width of the mass dis- 



tribution is consistent with the single star evolutionary 
scenario or the Mengel et al. binary evolution scenario 
(No. II in the introduction) for the origin of EHB stars, 
but can hardly be explained within the framework of 
the binary scenario (I) and the merger scenario (III) 
listed in the introduction, since both scenarios would 
predict a broader mass distribution with a mean mass 
below 0.5 Mq. 

The low masses as well as the low surface gravities we 
find for the BHB stars (brighter than the gap) remain 
a puzzle. Lower than expected gravities for BHB stars 
were also found recently by Dixon et al. ( |1996 ) in the 
globular cluster NGC 1904. Figure 10 illustrates this 
phenomenon with data for the globular clusters M 5, 
M 15, M 92, NGC 6397 and NGC 288 in addition to 
the data of this paper. As can be seen there the mass 
distribution of the BHB stars differs considerably from 
those of the sdB stars analysed here. The fact that the 
mean masses of BHB stars in globular clusters are sig- 
nificantly and systematically too low has already been 
noted in papers I and II. Most of the possible expla- 
nations of this finding (see papers I and II), however, 
are ruled out by the mass distribution of the sdB stars 
presented here: Non-solar helium abundances, peculiar 
metallicities due to diffusion, helium stratification, sys- 
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tematic errors in logg and/or distance moduli should 
affect the sdB stars in the same way as the BHB stars, 
(iv) Further insight into the nature of BHB and EHB stars 
may come from their radial distribution. We therefore 
investigated the radial distribution of all stars in Buo- 
nanno et al. (1986) that are fainter than V = 13™4 
and bluer than B— V = — 0™2. This sample, consisting 
of 214 stars, was divided into two groups: BHB stars 
with V < 16 (142), and EHB stars with V > 16 (70), 
which combines all sdB stars. We then counted the 
stars in radial bins of 50" and normalized the numbers 
to the total numbers of stars in each group as well as 
to the area covered by each bin. The resulting radial 
distributions show no significant difference, supporting 
a similar evolutionary history and mass distribution, 
in contrast to the situation described above. 

While the idea that there is something wrong with the 
model atmospheres, leading to too low values of logg for 
the temperature range between 10000 K and 20000 K 
sounds rather intriguing and would also explain the loci 
of the BHB stars in the log Tcff-log g-diagram (cf. Fig. 9 
of Paper I), it could not explain the findings of Paper II, 
where the physical parameters of the stars fit very well 
to theoretical expectations, but result nevertheless in too 
low masses. The remaining possibility, i.e. that the BHB 
masses are indeed that low puts severe problems to stellar 
evolutionary models, since up to now no scenario is known 
to produce stars in that temperature range with masses 
of about 0.3 M0. 
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Fig. 10. The resulting masses plotted in Fig. 9 compared to 
masses of BHB stars in other cluster. The BHB data are taken 
from Paper I (M 15), II (NGC 6397), and Crocker et al. (1988, 
M 5, M 92, and NGC 
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